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The 3’ end of the 20-kb genome of the Mebus strain of bovine enteric coronavirus (BCV) was copied into cDNA and 
cloned into the Pstl site of the pUC9 vector. Four clones from the 3’ end of the genome were sequenced either 
completely or in part to determine the sequence of the first 2451 bases. Within this sequence were identified, in order, 
a 3’-noncoding region of 291 bases, the gene for a 448-amino acid nucleocapsid protein (N) having a molecular weight 
of 49,379, and the gene for a 230-amino acid matrix protein (M) having a molecular weight of 26,376. A third large 
open reading frame is contained entirely within the N gene sequence but is positioned in a different reading frame; it 
potentially encodes a polypeptide of 207 amino acids having a molecular weight of 23,057. A higher degree of amino 
acid sequence homology was found between the M proteins of BCV and MHV (87%) than between the N proteins 
(70%). For the M proteins of BCV and MHV, notable differences were found at the amino terminus, the most probable 
site of O-glycosylation, where the sequence is N-Met-Ser-Ser-Val-Thr-Thr for BCV and N-Met-Ser-Ser-Thr-Thr for 
MHV. BCV apparently uses two of its six potential O-glycesylation sites. © 1987 Academic Press, Inc. 


INTRODUCTION 


The bovine enteric coronavirus (BCV) is one cause 
of severe enteritis in calves and may be responsible 
for as much as one-quarter of all deaths due to this 
disease (House, 1978). Vaccines produced from cell 
culture-attenuated strains of virus have failed to be 
completely protective. Before attempting to develop 
vaccines by recombinant DNA that may have improved 
usefulness, it is imperative that the genes and gene 
products responsible for inducing protective immunity 
be thoroughly characterized. Toward this end, and for 
the purpose of determining the function of individual 
proteins in coronavirus replication, we have begun to 
clone and sequence the BCV genome. 

BCV is known to possess a single-stranded, non- 
segmented, polyadenylated RNA genome of approxi- 
mately 20 kb (Guy and Brian, 1979; Lapps and Brian, 
1985). The total number of genes encoded by the ge- 
nome is not known, but presumably, because of its 
close antigenic relatedness to the well-characterized 
mouse hepatitis coronavirus, BCV will be similar to 
MHV in the number and arrangement of genes on its 
genome. One striking dissimilarity between BCV and 
MHV, however, is the possession by BCV of a fourth 
major structural protein, the 140-kDa hemagglutinin 
protein that comprises two disulfide-linked subunits of 
65 kDa (Hogue et a/., 1984; King and Brian, 1982; King 
et al., 1985). Questions therefore arise concerning not 
only its origin, function, and role in inducing protective 
immunity, but also the location of the hemagglutinin 


'To whom requests for reprints should be addressed. 


47 


gene on the genome and the resulting divergence from 
the MHV genome structure. 

In this paper we describe experiments that begin to 
examine the BCV genome by cDNA cloning and DNA 
sequencing. Within the 3’ 2451-base sequence we find 
a gene map that parallels that for MHV. We report the 
primary structure for the N and M genes and their de- 
duced amino acid sequences. Structural comparisons 
with other coronavirus N and M sequences are made 
and some conserved structural domains are identified. 


MATERIALS AND METHODS 


Virus and cells 


The Mebus strain of bovine coronavirus (BCV) was 
plaque purified and grown on the human rectal tumor 
(HRT) cell line as previously described (Hogue et a/., 
1984; Lapps and Brian, 1985). 


Radiolabeling of viral proteins and 
purification of virus 


Confluent monolayers of cells grown in 150-cm? 
flasks were infected with a multiplicity of approximately 
1 PFU per cell. After 1.5 hr adsorption at 37°, inoculum 
was removed and 15 ml of the appropriate medium 
and radioisotope was added. Viral polypeptides were 
labeled by adding 400 uCi 2H-labeled essential amino 
acids (150-200 mCi/mg; ICN) per flask in medium con- 
taining 10% normal essential amino acid concentration 
and 2% fetal calf serum (Sterile Systems, Inc.). Viral 
glycoproteins were labeled by adding 400 uzCi of 
[PH] glucosamine (5-15 Ci/mmol, ICN) per flask to me- 
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dium containing 5% fetal calf serum. Virus was har- 
vested and purified by isopycnic sedimentation in con- 
tinuous Sucrose gradients as previously described 
(Hogue et a/., 1984; Lapps and Brian, 1985). 


Polyacrylamide gel electrophoresis 
and immunoblotting 


The discontinuous buffer gel system of Laemmli 
(1970) was used as previously described (Hogue et a/., 
1984). For examining intracellular proteins, whole cell 
lysates were prepared by sonication. Cells in 60-mm 
petri dishes were washed twice with cold phosphate- 
buffered saline (PBS), scraped into cold PBS, and pel- 
leted by centrifugation at 2000 rpm. The cell pellet was 
suspended in 100 yl sterile distilled water, sonicated 
for 10 sec in a bath sonicator, and stored at —80°. For 
inhibitor studies, tunicamycin (Sigma) was used at a 
final concentration of 1.2 or 12 nM and monensin (Cal- 
biochem) was used at a final concentration of 1.0 uM. 
Tunicamycin or monensin was added to cells imme- 
diately after virus adsorption and was incubated with 
the ceils for a total of 24 hr, the time of cell lysate 
preparation. For electrophoresis, equal volumes of cell 
lysate and double-strength sample treatment buffer 
were mixed and heated at 100° for 5 min prior to elec- 
trophoresis. Unit strength sample treatment buffer is 
0.125 M Tris—HCl (pH 6.8)-4% sodium dodecyl sulfate— 
5 M urea. For immunoblotting, a modified method of 
Towbin et a/. (1979) was used as previously described 
(Hogue et a/., 1984). The preparation of rabbit anti- 
serum against individual BCV proteins was previously 
described (Hogue et a/., 1984). Monoclonal antiserum 
to human coronavirus 0C43 M protein, which also rec- 
ognizes BCV M protein, was a gift from J. Fleming, 
University of Southern California. 


Purification of genomic RNA 


Virus was purified from clarified supernatant fluids 
as described above. One-tenth of the virus preparation 
was labeled with [PH]uridine (400 Ci/mmol, ICN), 20 
uCi/ml, in order to follow RNA purification. Viral RNA 
was extracted using the proteinase K-SDS method 
(Lapps and Brian, 1985) and pheno!—chloroform-iso- 
amyl alcohol extraction and was ethanol precipitated 
after adding sodium acetate. Because subgenomic 
RNA species are incorporated into BCV virions (Lapps 
and Brian, 1985), full-length genomic RNA to be used 
for cDNA cloning and for making probe for colony 
screening was selected by rate-zonal sedimentation 
on preformed 5-mi linear gradients of 15 to 30% su- 
crose (w/w) made up in TNE (10 mM Tris-HCl (pH 7.5), 
100 mM NaCl, 1 mM EDTA)-0.1% SDS. RNA was dis- 
solved in water and sedimented for 1.5 hr at 110,000 
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g, 25°. RNA sedimenting faster than mammalian 28 S 
ribosomal RNA was recovered by ethanol precipitation. 


cDNA cloning of the 3’ end of the BCV genome 


BCV genomic RNA was cloned using a modified 
method of Gubler and Hoffman (1983). First-strand 
synthesis was carried out in a volume of 50 yl contain- 
ing 50 mM Tris-HCl (pH 8.1), 148 mM KCI, 8 mV 
MgCl., 1 mM DTT, 2 mM each of the four dNTPs, 10 
uCi [8*P]dCTP (3000 Ci/mmol, ICN), 15 units RNasin 
(Promega), 50 pmol oligo(dT 12-18), 3 wg BCV RNA, 20 
U reverse transcriptase (Seikagaku) for 1 hr at 37°, and 
the reaction was stopped by adding 2.5 ul 0.5 M EDTA. 
BCV RNA was heated to 100° for 5 min and quick 
cooled to 37° immediately before its addition to the 
reaction. Reaction products were extracted with phe- 
nol—chloroform—isoamy! alcohol and ethanol precipi- 
tated after adding ammonium acetate. 

Second-strand synthesis was carried out as de- 
scribed by Gubler and Hoffman in 100 ul containing 20 
mM Tris—HCl (pH 7.5), 5 mM MgCle, 10 mM (NH,)2S04, 
10 mM KCl, 0.15 MMB-NAD, 40 pM dNTPs, 8.5 U/ml 
Escherichia coli RNase H, 230 U/ml DNA polymerase 
1, 10 U/mi DNA ligase, and all of the product from first- 
strand reaction. Free nucleotides were removed by 
three cycles of ethanol precipitation of the reaction 
product and the total quantity of product was estimated 
from the amount of radiolabeled first strand that re- 
mained. 

Double-stranded cDNA was homopolymer tailed es- 
sentially by the method of Roychoudhury and Wu 
(1980). The following were added to the dried DNA in 
order: 20 wl 10 cacodylate—Tris buffer (1.4 4 K-cac- 
adylate, 0.3 M Tris-HCI (pH 7.6)), 4 nl 5 mM DTT, 3 ul 
10 mM dCTP, 162 pl HO, 2 wl 100 mM CoCl., 50 
uCila — **P]dCTP (+3000 Ci/mol, ICN) in 5 ul, and 16 
units of TdT in 2 yl. The reaction was carried out at 
37° for 1 min and then stopped by adding 10 wl 0.5 VM 
EDTA and 2 ul 10% SDS, and the product was phenol- 
chloroform—isoamy! alcohol extracted and ethanol pre- 
cipitated. This reaction assumed an average size of 1 
kbp for the ds cDNA and was designed to give an av- 
erage of 15 dCMP residues per 3’ end of dsDNA. 

C-tailed ds cDNA was annealed to G-tailed, Psti-lin- 
earized pUC9 vector (PL Biochemicals) and E. coli strain 
JM103 was transformed by the method of Hanahan 
(1983) using a total concentration of DNA of less than 
0.1 pg/ml. 


Identification of large inserts containing 3’-specitic 
BCV sequences 


Cells containing recombinant plasmids were ob- 
served as white colonies on YT agar plates that con- 
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tained 100 wg ampicillin/ml, 1 mM IPTG, 0.004% X-gal 
and were transferred to nitrocellulose (Millipore, HATF) 
and probed with random-primed cDNA copied from 
BCV genomic RNA (Maniatis et a/., 1982). **P-labeled, 
random-primed cDNA was synthesized as described 
above for the oligo(dT)-primed reaction except that 
dNTP concentrations were 2.5 uM each, 0.2 wg RNA 
was used, and oligo(dT) was replaced by 20 ug frag- 
mented calf thymus DNA. Colonies yielding strong sig- 
nals were analyzed for plasmid size and inserts of 1.0 
to 4.1 kb (the largest) were further analyzed by Southern 
hybridization with %P-labeled poly(dT) to detect 
poly(dA). *?P-labeled poly(dT) probe was prepared as 
described above for the oligo(dT)-primed reaction ex- 
cept that 50 pmol oligo(dT) - poly(rA) (PL Biochemicals) 
replaced the RNA. Alkali-treated [°?P]poly(dT) probe 
was incubated for hybridization at 37° for 12 hr, then 
at 20° for 36 hr, and blots were washed in 2x SSC, 
0.1% SDS at 20°. 


DNA sequencing and sequence analysis 


Plasmids were purified by alkaline lysis and cesium 
chloride centrifugation as described by Maniatis et a/. 
(1982), and restriction endonuclease mapping was 
done as described by Smith and Birnstiel (1976) using 
plasmids that were labeled at the Sa/l site within the 
multiple cloning linker region. Restriction fragments 
end-labeled with °*P were isolated and sequenced by 
the method of Maxam and Gilbert (1980). Many end- 
labeled fragments of less than 700 bases were first 
strand-separated before sequencing (James and Brad- 
shaw, 1984). Sequences were analyzed with the aid 
of the program developed by Queen and Korn (1984) 
marketed as part of the Beckman Microgenie program, 
March 1986 version (Beckman Instruments, Inc.). 


RESULTS 


cDNA cloning and sequencing of four clones from 
the 3’ end of the genome 


Starting material for cDNA cloning was approximately 
3 wg of rate-zonally purified genomic RNA obtained from 
500 ml of tissue culture medium. An estimated 70 ng 
of ds cDNA was generated and from this 670 white 
colonies were obtained. By colony screening, 89 col- 
onies gave a strong signal to [*7P]cDNA prepared from 
genomic RNA, and of these, 9 had inserts ranging from 
1.2 to 4.1 ko as determined by agarose gel electro- 
phoresis of linearized plasmids. The 9 clones were fur- 
ther analyzed to determine their restriction enzyme 
maps and poly(A) content. Only one of the clones, an 
insert of 1.2 kb identified as clone CB9, reacted strongly 
under hybridization conditions by Southern blotting to 
32P_labeled oligo(dT). Three other clones identified as 
MNS (2.1 kb}, MAS (2.8 kb), and MA7 (4.1 kb) were 


found to contain sequences that overlap with CB9 on 
the basis of hybridization and restriction endonuclease 
maps (data not shown). 

The orientation of all four clones in reference to the 
20-kb virus genome and the restriction enzyme sites 
used for sequencing are illustrated in Fig. 1. Our ori- 
entation presumes polyadenylation at only the 3’ end 
of the genome and this is based on the documented 
3’ polyadenylation site in the avian infectious bronchitis 
virus and mouse hepatitis virus genomes (Lai et al., 
1981; Stern and Kennedy, 1980). The strategy used 
for sequencing is described in the legend to Fig. 1. 
Initially, clone MN3 was sequenced completely and 
was found to contain all of the N and part of the M 
sequences. Greater than 98% of the sequence con- 
taining the complete N gene was determined either by 
sequencing both strands of clone MN3 or by repeated 
sequencing of the same strand using different methods 
of end labeling. Some of the sequences were confirmed 
from subclones of MN3. To complete the sequencing 
of the M gene, clone MA5 was sequenced from its 
second Ddel site and parts of MA7 were sequenced 
as described in Fig. 1. The total sequence of the M 
gene was determined by sequencing both strands of 
DNA and by repeated sequencing of some fragments 
using different methods of end labeling. 

The total nucleotide sequence of 2451 bases from 
the 3’ end of the genome and the deduced amino acid 
sequences for the three largest open reading frames 
contained in this sequence are illustrated in Fig. 2. 

All possible translation products were deduced for 
both virus-sense RNA and virus complementary-sense 
RNA (Fig. 3) because of the precedent that some single- 
strand RNA virus genomes are of ambisense polarity 
(Auperin et a/., 1984). RNA complementary to corona- 
virion-sense RNA could therefore theoretically function 
as mRNA. Figure 3 illustrates the three largest open 
reading frames found in the virus-sense sequence, 
each having the proper initiation codon and each being 
preceded by a termination codon. These are labeled 
N (for nucleocapsid protein), M (for matrix protein), and 
IORF (for ‘“‘internal’’ open reading frame). Within the 
virus complementary-sense RNA, two open reading 
frames of greater than 250 bases exist, beginning at 
approximate nucleotide positions 1220 and 1880 in 
the first reading frame. The significance of these is un- 
known. 

The largest open reading frame predicts a protein 
having properties of the nucleocapsid protein 

The largest open reading frame extends from base 
817 through base 2160 and predicts a 448-amino acid 
protein of 49,379 mol wt. We conclude this to be the 
coding sequence for the nucleocapsid protein (N) for 
the following reasons: (i) The only BCV protein de- 
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scribed to date that approaches this size is the 52-kDa virus strains A59 and JHM (Armstrong et a/., 1983; 
phosphorylated nucleocapsid protein (King and Brian, Skinner and Siddell, 1983). 


1982). (ii) The predicted protein is basic, a property The N gene for BCV shares other properties with the 
expected of nucleic acid-binding proteins. Fifty-nine N gene of MHV. (i) It is rich in serine. Forty-two residues 
(13%) of the amino acids are basic whereas 43 (10%) of serine make it the most abundant amino acid. (ii) It 
are acidic, giving the protein a net charge of +16 at is flanked on its 5’ side by the gene for the M protein, 
neutral pH. (iii) The amino acids encoded by this se- and it is flanked on its 3’ side by a noncoding region of 


quence share extensive (70%) sequence homology 291 bases, only 3 bases fewer than that for MHV JHM. 
with the N protein of the closely related mouse hepatitis (iii) The intergenic sequence between the M and N 
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Fic. 1. Sequencing strategy used to obtain BCV genomic sequences containing the N and M genes. (A) Strategy for obtaining the N gene 
sequence. Clone MN3 was sequenced completely and clones CB9 and MA5 were sequenced in part. The internal Ddel, Pstl, Sau3A |, and 
Xbal sites derived by restriction endonuclease mapping and the Hindlll and Sail sites in the multiple cloning region of the pUC9 vector were 
the sites used for DNA sequencing. @, | and @ indicate sites labeled at the 5’ end using polynucleotide kinase for clones MN3, CB9, and 
MAS, respectively. @, O, and @ indicate sites labeled at the 3’ end using reverse transcriptase, and the appropriate labeled deoxynucleotide 
triphosphate, for clones MN3, CB9, and MAS, respectively. indicates sites labeled at the 3’ end using radiolabeled cordycepin and terminal 
transferase for a subclone of MN3. (B) Strategy for obtaining the M gene sequence. Parts of clones MN3, MA5, and MA7 were sequenced 
beginning with the second Doel site from the 5’ end of clone MA5.i™, Z, and [indicate restriction sites that are 5’ end labeled with polynucleotide 
kinase for clones MN3, MA5, and MA7, respectively. @ and @ indicate sites that were 3’ end labeled using reverse transcriptase and the 
appropriate labeled deoxynucleotide triphosphate for clones MN3 and MAS, respectively. Uniquely labeled molecules for sequencing were 
obtained from gels after electrophoretic separation of restriction endonuclease-treated, end-labeled fragments or after strand separation. The 
orientation of the clones in the pUC9 vector are as shown except for CB9 which is inverted; i.e., the poly(A) end is next to the Hindill site of 
the vector. pUC9 sequences are indicated as a bold line at the end of the restriction endonuclease maps. 
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genes is very similar. Beginning with the first base fol- 
lowing the M gene termination codon, the sequence 
is UAUCUAAACUUUAAGG for BCV, and UCUAAAC- 
UUUAAGG for MHV. (iv) The consensus sequence sur- 
rounding the initiation codon for the N gene, AG- 
GAUGU, is the same, and is a preferred sequence for 
translation initiation among eukaryotic messenger 
RNAs (Kozak, 1983). 


The second largest open reading frame predicts a 
protein having properties of the matrix protein and 
also identifies potential O-glycosylation sites 


The second largest open reading frame extends from 
base 112 through base 804 and predicts a 230-amino 
acid protein having a molecular weight of 26,376 (Fig. 
2). This protein has extensive amino acid homology 
with the M protein of MHV A59 strain (Armstrong et 
al., 1984), as expected from its close antigenic relat- 
edness (Hogue et a/., 1984), and is therefore the ap- 
parent BCV counterpart (Fig. 4). By maximum alignment 
of the proteins, 200 of the amino acids (>86%) are the 
same as 200 of the 228 in the MHV sequence, and 
another 16 (7%) represent conservative changes. Be- 
cause of the strong similarity in structure between the 
BCV and MHV M proteins, the BCV M protein can be 
expected to have a similar topology with respect to the 
virion envelope (Rottier et a/., 1986). Namely, the central 
portion of the molecule would be expected to span the 
membrane three times, with approximately 28 amino 
acids (26 for MHV) at the amino terminus being external 
to the virion and approximately 100 amino acids at the 
carboxy terminus being internal to the virion. The pro- 
tein is slightly basic, having a net charge of +9 at neutral 
pH. The basic amino acids are clustered in the carboxy 
terminal 40% of the protein. Within the carboxy terminal 
100 amino acids are 14 of the 20 basic amino acids 
and 6 of the 14 acidic amino acids, giving this region 
of the molecule a net charge of +8. It is therefore rea- 
sonable to expect that this part of the molecule might 
be interacting with the negatively charged RNA as sug- 
gested (Sturman et a/., 1980) or possibly with an acidic 
portion of the N protein to contribute to a direct inter- 
action between the M and N molecules. We predict 
the latter occurs on the basis of a 1:1 molar ratio be- 
tween the M and N proteins in BCV (King and Brian, 
1982). One hundred thirteen (49%) of the amino acids 
are hydrophobic and the distribution of hydrophobic 
amino acids is nearly identical to that for the MHV M 
protein. 


Evidence for two O-linked oligosaccharides 
per M molecule 


The M proteins of BCV and MHV were together the 
first viral glycoproteins shown to possess O-linked oli- 


gosaccharides (Holmes et a/., 1981; Niemann and 
Klenk, 1981). The character of the oligosaccharides, 
however, has been described only for MHV A59 (Nie- 
mann et a/., 1984). Our data suggest there may be up 
to two O-linked oligosaccharides per BCV M molecule. 
First, three separate species of M (gp26) molecules 
were identified from purified BCV when radiolabeled 
proteins were resolved by electrophoresis (Fig. 5). 
These were also observed, but less clearly resoived, 
when identified by immunoblotting with M-specific 
polyvalent antiserum (Fig. 5) or with M-specific mono- 
clonal antiserum (data not shown). They have apparent 
molecular weights of 22K, 24K, and 26K and their ap- 
pearance is consistent with the notion that the 22K 
species is the unglycosylated precursor and 1 or 2 oli- 
gosaccharide chains, each contributing approximately 
2 kDa toward the molecular weight (Klenk and Rott, 
1981), are added to assemble a 24- and 26-kDa spe- 
cies, respectively. Second, only three species of M 
protein were resolved in lysates of infected cells by 
immunoblotting and neither their sizes nor relative 
amounts were altered by tunicamycin, an inhibitor of 
N-glycosylation but not O-glycosylation. Tunicamycin 
does, however, inhibit the glycosylation of gp190 (the 
peplomeric protein for which the virion-associated 
subunits are gp120 and gp100) and gp140 (the hem- 
agglutinn) (Hogue and Brian, in preparation). In the 
presence of radiolabeled glucosamine only the 24- and 
26-kDa species were labeled (Fig. 5). The fact that mo- 
nensin, an inhibitor of Golgi function and hence O-gly- 
cosylation, diminishes the amount of the 24- and 26- 
kDa species and enhances the relative abundance of 
the 22-kDa species strengthens the notion the M gly- 
cosylation is O-linked (Niemann et a/., 1982). 

Assuming the BCV M protein is glycosylated in the 
region external to the virion envelope, i.e., within the 
first 28 amino acids of the N terminus, then the serine 
residues at positions 2 and 3 or the threonine residues 
at positions 5, 6, 12, and 14 are potential sites for O- 
glycosylation (Fig. 4). If, as presumed for MHV, the gly- 
cosylation sites are primarily within the N-terminal N- 
Met-Ser—Ser—Thr—Thr sequence, a region identical to 
the glycosylated region of glycophorin A, then the se- 
quence per se may not be an absolute requirement for 
glycosylation since the N terminal sequence for BCV 
is N-Met-Ser—Ser—Val-Thr-Thr. 

The discrepancy between the observed molecular 
weight of 22 kDa for the unglycosylated polypeptide 
and the molecular weight of 26,376 deduced from se- 
quence data could be explained by a strong tendency 
of the hydrophobic regions of the M protein to remain 
in close proximity, even in the presence of SDS, giving 
rise to more rapidly migrating globular molecules. Cer- 
tainly such behavior would explain the self-aggregation 
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. 30 60 90 120 

AGGACTGTCCCCTICTATTTATGTGTTTAATAGAGGTAGGCAGTTTTATGAGTITTACAACGATGTAAAACCACCAGTICTTGATGTGGATGACGTTTAGTTAATCCAAACATTATGAGT 
M S 

150 180 . . 210 240 


AGTGTAACTACACCAGCACCAGTTTACACCTGGACTGCTGATGAAGCTATTAAATTCCTAAAGGAATGGAACTTTTCTITGGGTATTATACTACTITTTATTACAATCATATTGCAATT? 
s v T T P A PV ¥ T WT AD E AT K F CK Ew N F SS LGIILLFIftItuLua9 F 


270 300 330 360 
GGATATACAAGTCGCAGTATGTTTGTTTATGTTATTAAGATGATCATTTTGTGGCTTATGTGGCCCCTTACTATCATCTTAACTATTTTCAATTGCGTGTATGCGTTGAATAATGTGTAT 
GY T S&S RS MFVY¥YvVdIK MI1ILWoLOMwWPL Tf ITtTobttitiFNne¢ecvyawunNwNvVvV Y 


390 420 450 480 
CTIGGCTITTCTATAGTTTTCACTATAGTGGCCATTATCATGTGGATIGTGTATTTTGTGAATAGTATCAGGTTGTTTATTAGAACTGGAAGTTGGTGGAGTTTCAACCCAGAAACAAAC 
L GF S$ IV FTI V ATIdiIM WwW IV ¥ F VN S IT ROU F I RTGS WW S F N P E TN 


510 . 540 570 600 
AACTTGATGTGTATAGATATGAAGGGAAGGATGTATGTTAGGCCGATAATTGAGGACTACCATACCCTTACGGTCACAATAATACGTGGTCATCTTTACATGCAAGGTATAAAACTAGGT 
NLM C IDM K G RM Y V R PITTED Y HT LTV TIT RGHOoLY MQGtI KL G 


630 660 690 720 
ACTGGCTATTCTTTGTCAGATTTGCCAGCTTATGTGACTGTTGCTAAGGICTCACACCTGCTCACGTATAAGCGTGGTTTTCTTGACAAGATAGGCGATACTAGTGGTTTTGCTGTTTAT 
TG YS LS DL PAYVTVA KV S HOUBULT YY K RGF LODK I GDTS GFAV ¥ 


750 . 780 810 840 
GTP'TAAGTCCAAAGTCGGTAATTACCGACTGCCATCAACCCAAAAGGGTTCTGGCATGGACACCGCATTGTTGAGAAATAATATCTAAACTTTAAGGATGTCTTTTACTCCTGGTAAGCAA 
v eK S K V GN ¥ RL PS TQK GS GM DTA LLR NN I M S F T P G K Q 

. 870 . 900 930 . 960 


TCCAGTAGTAGAGCGTCCTTTGGAAATCGTTCTGGTAATGGCATCCTTAAGTGGGCCGATCAGTCCGACCAATCTAGAAATGTTCAAACCAGGGGTAGAAGAGCTCAACCCAAGCAAACT 
5 S$ S$ RAS FGNRSGNGtHILKWNAD QS DQSE RN V QTR GR RA QP K OQ T 
MAS LS GPI SS PT N LEM F K P GV E E LN PS K OL 


990 1020 . 1050 . 1080 
GCTACTTCTCAGCTACCATCAGGAGGGAATGTTGTACCCTACTATTCTTGGTTCTCTGGAATTACTCAGTTTCAAAAAGGAAAGGAGTTTGAATTTGCAGAGGGACAAGGTCTGCCTATT 
ATS QLPS GGNVV PY ¥ S WF S GIT Q FQ K G K BE F E F A EGQ GV PI 

LLLSs YY HQ@EGM LY PTI LG S&S L ELLs F K K ER S LNLQRoODK Vc¢cL UL 


1110 . 1140 . 1170 1200 
GCACCAGGAGTCCCAGCTACTGAAGCTAAGGGGTACTGGTACAGACACAACAGACGTTCTTTTAAAACAGCCGATGGCAACCAGCGTCAACTGCTGCCACGATGGTATTTTTACTATCTT 
A P GV PAT EA K GY WY RHNRRS FK TAD GN QR QL LP RW Y¥ F Y¥ ¥ L 

HO -E S Oo D. i K DR GT G6 f DT T DV LL LK QO P.M A TS VON eC. CO Be DooG TUR Sf aT. Gh 


. 1230 . 1260 . 1290 1320 
GGAACAGGACCGCATGCCAAAGACCAGTATGGCACCGATATTGACGGAGTCTTCTGGGTCGCTAGTAACCAGGCTGATGTCAATACCCCGGCTGACATTCTCGATCGGGACCCAAGTAGC 
GTGPHAK DQYG6G?DIDG YF WV AS N Q AD VN T PA OD Tbe Do Re UD. (Be S48 

EOobDR*M PK TS MAPILT ES S GCGSLEVTRLUM s I P RLUTF SI G6TQV A 


. 1350 . 1380 1410 . 1440 
GATGAGGCTATTCCGACTAGGTTTCCGCCTGGCACGGTACTCCCTCAGGOTTACTATATTGAAGGCTCAGGAAGGTCTGCTCCTAATTCCAGATCTACTTCACGCGCATCCAGTAGAGCC 
DEA IPTR*FPPGTVL POG Y ¥ IEG S GR S APN S_R $s T S RAS S RA 

MROLFROGFROLaARYS LRVTIL kK AQ EGLLL I P D L LH AH P V €E P 


1470 . 1500 1530 1560 
TCTAGTGCAGGATCGCGTAGTAGAGCCAATTCTGGCAACAGAACCCCTACCTCTGGTGTAACACCTGATATGGCTGATCAAATTGCTAGTCTTGTTCTGGCAAAACTIGGCAAGGATGCC 
SSAGS RS RANSGNRTPTSGCOVTPDOMADQIAS LEVLAK LG KDA 

L VQ oD RV V E P ILA T E PL PE Vv 


. 1590 . 1620 . 1650 . 1680 
ACTAAGCCACAGCAAGTAACTAAGCAGACTGCCAAAGAAATCAGACAGAAAATTTTGAATAAGCCCCGCCAGAAGAGGAGCCCCAATAAACAATGCACTGTTCAGCAGTGTTTTGGGAAG 
TK POove?xK OT AK ET RQK IL NK PROK RS PNK Q ET V QQ CF GK 


* 1710 . 1740 ‘ 1770 . 1800 
AGAGGCCCCAATCAGAATTTTGGTGGTGGAGAAATGTTAAAACTIGGAACTAGTGACCCACAGTICCCCATTCTTGCAGAACTCGCACCCACAGCTGGTGCGTTTTTCTTTGGATCAAGA 
R GPNQN*FG6GeGeEwm”MULkKiOULGtTt s DPQ EF PIL A E L AP T AGA F F F GS R 


1830 1860 1890 1920 
TTAGAGTTGGCCAAAGTGCAGAATTTGTCTGGGAATCTTGACGAGCCCCAGAAGGATGTTTATGAATTGCGCTATAATGGTGCAATTAGATTTGACAGTACACTTTCAGGTTTTGAGACC 
L ELA KVYQNULUS GNLOD EP QK DV Y EL R ¥Y N GA IRF ODS TLS GF ET 


. 1950 S 1980 . 2010 2040 
ATAATGAAGGTGTTGAATGAGAATTTGAATGCATATCAACAACAAGATGGTATGATGAATATGAGTCCAAAACCACAGCGTCAGCGTGGTCAGAAGAATGGACAAGGAGAAAATGATAAT 
I MK VLNENLNA ¥ Q@=QQD GM MNM S PK PQ@gRQRGOQKNGQ GEN ODN 


2070 . 2100 . 2130 2160 
ATAAGTGTTGCAGCGCCTAAAAGCCGTGTGCAGCAAAATAAGAGTAGAGAGTTGACTGCAGAGGACATCAGCCTTCTTAAGAAGATGGATGAGCCCTATACTGAAGACACCTCAGAAATA 
Is Vv AA PK S RVQOQNK S$ R ELTA EDI: S LLK K MD EP ¥ TE DT S E I. 


2190 2220 2250 2280 
TAAGAGAATGAACCTTATGTCGGCACCTGGTGGTAAGCCCTCGCAGGAAAGTCGGGATAAGGCACTCTCTATCAGAATGGATGTCTIGCTGCTATAATAGATAGAGAAGGTTATAGCAGA 


2310 2340 2370 2400 
CTATAGATTAATTAGTTGAAAGTTTTGTGTGGTAATGTATAGTGTTGGAGAAAGTGAAAGACTTGCGGAAGTAATTGCCGACAAGTGCCCAAGGGAAGAGCCAGCATGTTAAGTTACCAC 


2430 
CCAGTAATTAGTAAATGAATGAAGTTAATTATGGCCAATTGGAAGAATCAC 
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Fic. 3. Schematic diagram of open reading frames obtained when translating the 3’ 245 1-base sequence of the BCV genome as either virus- 
sense RNA or virus complementary-sense RNA. Vertical bars above the line represent the first methionine codon that could serve as the 
initiation site for translation. Vertical bars below the line represent termination codons. M: sequence of the matrix protein gene. N: sequence 
of the nucleocapsid protein gene. IORF: sequence of an internal open reading frame within the N gene. 5’ to 3’ orientation is always left to 


right. 


of the M protein observed so frequently (Fig. 5; Hogue 
et al., 1984). 


A large open reading frame found internal to the N 
gene, but in a different reading frame, predicts a 
207-amino acid protein having a molecular 

weight of 23,057 


The nucleotide sequence from base 878 through 
base 1498 in the second reading frame encodes a 207- 
amino acid protein of 23,057 (Fig. 2). This protein is 
hypothetical since we have no proof yet of its existence. 
The protein has a net charge of +1 at neutral pH and 
is moderately hydrophobic since 79 (38%) of its amino 
acids are hydrophobic. The hydrophobic amino acids 
are spread somewhat evenly throughout the protein 
except at the carboxyl terminus, where there are 
enough to make this part of the protein a potential 
membrane anchor region. Twenty-seven of the terminal 
41 amino acids (66%) are hydrophobic. The existence 
of the protein cannot be ruled out on the basis of the 
consensus sequence (GUAAUGG) surrounding its ini- 
tiation codon since it is one commonly used, being 
found at the initiation site of 18% of all eukaryotic 
mRNAs catalogued (Kozak, 1983), nor can it be ruled 


out on the basis of codon usage since it is similar to 
that used for the N and M proteins. 


DISCUSSION 


We present the first nucleotide sequence data avail- 
able for BCV or for any member of the hemagglutinating 
mammalian coronavirus subgroup which includes the 
human respiratory coronavirus OC43 and the porcine 
hemagglutinating encephalitis virus. Despite the fact 
that BCV has the hemagglutinin structural protein that 
is missing on MHV A59 (Hogue et a/., 1984; King et 
al., 1985), it shares membership with MHV in one of 
the four major antigenic subgroups of coronaviruses 
(Pedersen et a/., 1978). Both the gene map and the 
primary sequence for that part of the BCV genome de- 
scribed in this paper reflect a close relatedness to MHV, 
consistent with patterns of shared antigenicity between 
the two viruses (Hogue et a/., 1984). Genome sequence 
divergence with regard to the hemagglutinin gene must 
therefore lie 5’-ward of this sequence. Both gene ar- 
rangement and primary sequence at the 3’ end of the 
genome, however, suggest a greater degree of evo- 
lutionary divergence from both the porcine transmis- 
sible gastroenteritis virus (TGEV) and avian infectious 


Fig. 2. The primary nucleotide sequence of the 3’ 2451 bases of the BCV genome and the deduced amino acid sequences for the three 
largest open reading frames, M, N, and IORF. The M open reading frame extends from base positions 115 through 804, the N open reading 


frame from 817 through 2160, and the “‘internal’’ open reading frame from 878 through 1498. The conserved intergenic sequences, cf Aaac, 
are underlined with a solid line. A 10-base sequence highly conserved within the 3’-noncoding region among coronaviruses (Kapke and Brian, 


1986) is identified with a double underline. 
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Fic. 4. Comparison of the deduced amino acid sequences for the M proteins of the bovine coronavirus (top) and the mouse hepatitis virus 
A539 strain (bottom). In the MHV sequence, amino acids matching with those of the BCV sequence are indicated by a hyphen; unmatched 
amino acids are named. Gaps were introduced to maximize alignment. Potential O-glycosylation sites on the amino terminus of the BCV M 
proteins are indicated by solid circles. 


bronchitis virus (IBV) than from MHV. TGEV has an open (Kapke and Brian, 1986); IBV has two open reading 
reading frame for a potential 9.1K protein positioned frames for proteins of 7.5K and 9.5K positioned be- 
between the 3’-noncoding sequence and the N gene tween the N and M genes (Boursnell and Brown, 1984). 


ANTI gp26 3H-aa ANTIgp26 __3H-gluN 
VIRION LYSATE VIRION LYSATE VIRION 
™M - + + a ee ee ee ee 
MON - - - - 
2-ME - - = = 


-26 
-24 


1 2 3 4 5 6 7 8 9 10 


Fic. 5. Identification of different glycosylation states of the M protein. Lane 1: Purified virus immunoblotted with anti-go26 serum. Lanes 2, 
3, and 4: Lysates of infected cells grown respectively in 12, 1.2, or 0 4M tunicamycin (TM), electrophoresed on the same gel as lane 1, and 
immunoblotted with anti-gp26 serum. Lanes 5 and 6: Purified virus radiolabeled with *H-labeled amino acids. Lanes 7 and 8: Lysates of infected 
cells grown in 0 or 1 #M monensin (MON) and immunoblotted with anti-go26 serum. Lanes 9 and 10: Purified virus radiolabeled with 
[H]glucosamine. Samples treated with 2-mercaptoethanol are indicated (2-ME). The polyacrylamide concentrations are 9% for lanes 1 through 
4, and 8% for lanes 5 through 10. 
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The N protein of BCV shows an overall amino acid 
sequence homology of 70% with both MHV A59 and 
MHV JHM (72% at the nucleotide level) (Skinner and 
Siddell, 1983; Armstrong et a/., 1983) but only 29% 
(37% at the nucleotide level) with the N protein of TGEV 
(Kapke and Brian, 1986), and 29% (43% at the nucleo- 
tide level) with the N protein of IBV (Boursnell et a/., 
1985). The degree of homology between the N amino 
acid sequences of BCV and MHV is not evenly distrib- 
uted throughout the gene, however. There are regions 
of up to 16 amino acid stretches, for example, that 
show less than 30% homology. Conversely there are 
regions of up to 69-amino acid stretches showing 
greater than 90% homology. A region of high homology 
among MHV (beginning at amino acid 86), IBV (begin- 
ning at amino acid 53), and TGEV (beginning at amino 
acid 53), and extending for 68 positions (Kapke and 
Brian, 1986), is also found in BCV (beginning at amino 
acid 83). Within this region there is 79% perfect ho- 
mology between BCV and MHV. Such regions of con- 
servation suggest that there are evolutionary pressures 
for retention of a specific function associated with this 
sequence. Other regions having similar chemical prop- 
erties but little primary sequence homology also sug- 
gest conserved functional domains. These include 
clusters of serine residues and clusters of basic and 
acidic amino acids. Assuming all coronavirus N proteins 


BASIC AMINO SERINE 
ACID POSITIONS POSITIONS 


ACIDIC AMINO 
ACID POSITIONS 


100 200 


are phosphorylated at only serine residues, as in the 
N protein of MHV (Stohiman and Lai, 1979), then “hot 
spots" for potential phosphorylation become apparent 
when the N protein sequences are compared (Fig. 6). 
By aligning the N proteins of MHV, BCV, IBV, and TGEV 
with the first amino acid of the conserved 68 amino 
acid region, three clusters of 3-12 serine residues in 
common among all viruses become apparent at BCV 
amino acid positions 40-70, 180-225, and 300-350. 
The major serine cluster region is at amino acid posi- 
tions 180-225. Cluster groups of 5 to 26 basic amino 
acids can be seen within 50 residues of the amino 
terminus, within the 68-amino-acid conserved region, 
between amino acid positions 200 and 300, and in a 
region extending between amino acids 50 and 25 from 
the carboxy terminus, but not within six positions of 
the carboxy terminus (Fig. 6). Clustering of acidic amino 
acids is less striking but clusters of 10 to 12 are ob- 
served within the last 100 bases of the carboxy ter- 
minus (Fig. 6). Such regions may indicate sites for pro- 
tein—nucleic acid or protein-protein interactions. 

The high degree of amino acid sequence homology 
between the M proteins of BCV and MHV (86%) con- 
trasts with the lower degree (70%) between the N pro- 
teins. The contrast becomes even more striking when 
amino acids of conserved nature are included, making 
the homology 93 and 79%, respectively, for the M and 


300 400 


AMINO ACID POSITION 


Fic. 6. Positions of serine residues, basic amino acids, and acidic amino acids on the N protein of MHV (A59), BCV (Mebus), IBV (Beaudette), 
and TGEV (Purdue), depicted in a way to show clustering. The amino terminus is positioned on the left and the proteins (depicted with no gaps 
in their sequence) are aligned by the first amino acid in the 68-amino-acid region of high homology which is indicated by the stippled line. 
Amino acid positions are indicated by vertical bars. Proteins are arranged in decreasing order of length. 
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N proteins. This contrast indicates either that structural 
constraints on the M protein are more rigid, resulting 
in a more limited evolution of this protein, or that there 
is a form of genetic exchange that has taken place 
between the two viruses. The notion that the M protein 
may be structurally constrained as a result of functional 
requirements is suggested by the conserved chemical 
features between the MHV and IBV M proteins in the 
absence of conserved primary structure. IBV is anti- 
genically unrelated to MHV, and the IBV M protein 
shares an amino acid sequence homology of only 35% 
(perfect match only using the same method of align- 
ment employed above) with that of MHV. Yet it shows 
an extremely similar hydrophobicity profile and thus an 
apparently similar membrane topology (Boursnell et a/., 
1984). That is, amino acid changes were conservative. 
The notion of genetic exchanges, similar to those ob- 
served for RNA viruses with segmented genomes, must 
be seriously considered in light of recent evidence that 
coronaviruses undergo high-frequency recombination 
(Makino et a/., 1986). The mechanism giving rise to 
coronavirus recombinants is unknown but may involve 
displacement of nascent RNA polymerase complexes 
from the negative-strand template of one parent with 
subsequent attachment to the negative strand of a 
second parent (Makino et a/., 1986). Recombination 
might therefore be expected between the closely re- 
lated BCV and MHV viruses if, by chance, they should 
replicate simultaneously in the same host. This most 
certainly would be expected if polymerase binding dur- 
ing the recombinational event involves the conserved 
intergenic sequences used to identify initiation sites for 
transcription (Baric et a/., 1985; Budzilowicz et a/., 1985; 
Makino et a/., 1986). 
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